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We present a search for supersymmetry (SUSY) in the fully hadronic final state with the CMS detector at 
the LHC. This final state contains at least two jets and a significant transverse energy imbalance due to the 
neutralinos escaping detection. The background to the all-hadronic signature arise from QCD multi-jet pro- 
duction, ttbar and electroweak boson+jet production. These backgrounds can be estimated by utilizing a 
data-driven approach thus enabling a possible discovery of SUSY in the early physics data. In addition to the 
generic search, we also present an analysis which utilizes a new quantity, bj, constructed exclusively from the 
transverse energies of the jets to effectively eliminate the QCD background. 



1. Introduction 

The supersymmetry [lj](SUSY) is a potential sym- 
metry between the fermions and the bosons of the 
standard model (SM) of particles proposing that each 
fermion has a bosonic superpartner and each boson 
has a fermionic superpartner. The main attraction 
of the SUSY models is that they solve the hierarchy 
and the fine-tuning problems of SM. Given that no 
sparticles have been observed it must be that SUSY 
is broken such that the sparticles are heavier than the 
particles. 

One of the most studied SUSY models is the min- 
imally supersymmetric standard model (MSSM) of 
particles. This model is an extension of SM postu- 
lating the conservation of the R-parity. This conser- 
vation law results in the pair prodcution of the sparti- 
cles and their decay into an odd number of sparticles. 
This means further that the lightest supersymmetric 
particle (LSP) is stable and interacts very weakly with 
the matter making it a dark matter candidate. 

Another postulate of MSSM is the soft SUSY break- 
ing. There are several theoretical approaches to this, 
but this report will mention only a gravity mediated 
SUSY breaking (mSUGRA) case. 

The strategy used in every search for SUSY is 
geared towards selecting events with very energetic 
observables: leptons, photons, jets, missing transverse 
energy (Et ) due to neutrinos and LSP. The event se- 
lection and the methods used in these searches Q are 
optimized on the simulation of a set of benchmark 
points defined by certain values of the SUSY parame- 
ters. These points are chosen such that a good cover- 
age of the phase space is obtained beyond the current 
experimental and theoretical limits. 



2. The CMS detector and event clean-up 

The SUSY searches described in this report are us- 
ing simulated 10 TeV proton-proton collisions followed 
by the simulation of the detector response to the re- 
sulting particles. The detector used in these studies 
is the general purpose CMS [|[ detector installed at 



LHC. The simulation of this detector is made using 
the software package GEANT4 Q that takes into ac- 
count a detailed description of the geometry of the 
detectors, the materials and the magnetic field. 

In order to perform any kind of analysis with the 
data obtained from this detector, all the observables 
(i.e. jets, electrons, muons, photons, Et) have to be 
well understood. Given its discriminating power in 
any search for SUSY, If^ needs special attention. Be- 
sides potential undetected new particles, the sources 
for are many: cracks in the detector, cosmic rays, 
beam halo, electronic noise, and even SM processes 
(e.g. semileptonic decay of hadrons). At CMS it 
has been studied that the effect of the cosmic rays 
and of the beam halo can be minimized by selecting 
events with significant average electromagnetic frac- 
tion [EEMF > 0.1) and average charged fraction 
(ECHF > 0.175). 
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3. Search for SUSY program at CMS 

The strategy employed by the CMS collaboration to 
search for SUSY follows three main directives: model 
independent classification of the searches based on the 
topology of the final state, data-driven estimation of 
the SM backgrounds, and optimization of the event 
selection based on a set of benchmark signal points. 

Depending on the physics objects content of the 
final state, there are eight main SUSY searches at 
CMS: exclusive n-jet, inclusive 3-jets, photon+^T, 
single lepton, same-sign dilepton, opposite-sign dilep- 
ton, trilepton, and dilepton+photon analysis. In this 
report, only the hadronic searches (i.e. exclusive n- 
jet and inclusive 3-jets) are discussed. For each of 
these analyses the estimation of the various SM back- 
grounds is performed using techniques that rely as 
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little as possible on the detector simulation. Some of 
these techniques are described briefly in Scction [3.1.11 
In order to enhance the sensitivity of these searches 
the selection criteria used in these searches are lightly 
optimized using a set of benchmark SUSY signal 
points not excluded by theoretical and other hadron 
colliders experimental constraints. The distribution 
of some of these points is illustrated in Fig. [1] for the 
mSUGRA scenario. In the following sections the de- 
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Figure 1: The mSUGRA benchmark points used by the 
CMS collaboration to optimize the selection criteria ap- 
plied in the SUSY searches. 



tails of the inclusive hadronic search and of the exclu- 
sive multi-jet search are discussed. 

3.1. Inclusive hadronic search 

In this analysis we look for an excess over the SM 
background in a final state that contains at least three 
jets and significant ■ The trigger used to collect 
the data contains requirements on the scalar sum of 
transverse energies (Et) of the jets (HT) and on the 
magnitude of the vectorial sum of the transverse mo- 
mentum (p T ) of the jets (MHT). 
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At Level 1 we ask for HT > 200 GeV while at the 
High Level Trigger (HLT) the HT > 350 GeV and 



MHT > 100 GeV. In addition to the event clean-up 
described in Section [2j the offline analysis requires at 
least three jets with Et> 30 GeV and pseudo-rapidity 
(rf) < 3 in absolute value. The electron or photon 
contamination of these jets is reduced by requiring 
the electromagnetic fraction of the jet energy to be < 
0.9. 

The QCD background can be significantly reduced 
by selecting a certain region in the S(f>i-5(f>2 plane de- 
fined such that i?i > 0.5 and i? 2 > 0.5, where 



#1(2) = v/^l(l) + (?T - <%(2)) 2 



(5) 



and 84>i(2) is the difference in the polar angles of the 
leading(subleading) jet and I^t ■ The contribution 
from top quark pair production as well as from W 
boson production is reduced by requiring no isolated 
tracks in the event. 

To reduce the overall SM background, the leading 
jet is required to have Et> 180 GeV and \r]\ < 1.7, 
while the second leading jet should have Et> 110 
GeV. For further background reduction we require Ifi^ 
> 200 GeV and M eff > 500 GeV, where M eff is the 
sum of tfi'Y and the EtoI the second, third and fourth 
jets. This is a prototypical event selection in this anal- 
ysis that has been optimized using MC samples gener- 
ated at 14 TeV collisions and described in 0. For less 
energetic collisions, the values of the kinematic cuts 
should be lower than shown above and are currently 
under study at CMS. 

3.1.1. SM background 

An important part of the searches for SUSY is the 
evaluation of the SM background. The SM processes 
have a much larger cross-section than the potential 
SUSY signal. For the LM1 benchmark point (mo = 60 
GeV, mi/a = 250 GeV, tan/? = 10, A = 0, /z > 0) the 
leading order (LO) cross-section at 14 TeV is about 
42 pb. The corresponding LO cross-sections for the 
most important SM backgrounds are: «5.6E10 pb for 
QCD events, ssl5E3 pb for events with a Z-boson 
produced in association with jets, and «800 pb for 
tt events. In order to minimize the systematic uncer- 
tainties due to potential differences between data and 
simulation, the strategy behind the evaluation of these 
backgrounds relies on data driven methods. Also with 
these methods there is no need to calibrate the simu- 
lation to the data. 

The most important SM background is due to QCD 
processes. Using the cuts on the angle between the 
leading jets and (see Section [3~Tt . about 80% of 
the QCD events can be rejected, while keeping about 
90% of the SUSY events. The remaining QCD content 
can be evaluated directly from data via two proposed 
methods: smearing method and ABCD method. 

The smearing method, developed at ATLAS 0], 
relies on the parameterization of a response func- 
tion from multijet events with high values and I^t 
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aligned with one of the jets. The response function is 
defined for each jet as 



R 
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and it used to smear the jets from events with low If^ 
. The smearing of the jets will result in artificially 
created ^rused to estimate the real distribution. 
The normalization is obtained from the multijet data 
events with low $t • 

The ABCD method, pursued at CMS, uses two mi- 
correlated observables (I^t variable and the minimum 
azimuthal angle between its direction and that of the 
three leading jets). The plane formed by these two 
variables is split in four regions: A, B, C and D, 
such that the SUSY signal is contained in region C. 
Given that the two observables are uncorrelated, the 
number of QCD events in region C can be derived 
from the number of QCD events in the other regions: 
C = DxB/A. For this method to work it is important 
to have little SUSY content in region A, B and D. 

Another important SM background is produced by 
Z — ► vv +jets events. This background is irre- 
ducible and there are few proposed ways to mea- 
sure it. One method relies on the measurement of 
^xfrorn Z — > fifi -fjets events where the muons are 
ignored. The normalization is set by taking into ac- 
count the theoretical ratio of cross-sections between 
the Z — y vv and Z — > /i/i processes. Another method, 
developed at CMS p, relies on 7+jets events with the 
photon being ignored in the event. This method bene- 
fits from a larger data sample, while the disadvantage 
is due to the uncertainty on the normalization. The 
plot in Figure [2] shows the comparison between the 
distribution from Z —y vv +jets MC sample and 
the corresponding distributions using 7+jets events. 
The Z — y vv +jets events background can also be es- 
timated from leptonic W decays where the lepton is 
ignored. This method [|| also benefits from a larger 
statistics compared to Z — > /i/i +jets, but less than 
7+jets events. 

3.1 .2. Sensitivity reach 

In Fig. [3] it is shown the 5-sigma discovery limit in 
the mSUGRA mo-nix « plane for the final state with 
at least three jets and Jf^ . These curves are shown 
assuming 100 pb _1 of good data collected by CMS for 
14 TeV collisions (light color) and 10 TeV collisions 
(dark) . The limits obtained from the Tevatron exper- 
iments (solid, dark) and from LEP (dashed lines) are 
also shown. The curve for the 10 TeV collisions has 
been obtained using the same event selection used for 
the 14 TeV case. We conclude that for this final state 
there is a very good chance of discovering SUSY with 
several tens of pb _1 of understood data. 
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Figure 2: The comparison between the distribution 
from Z —y vv +jets MC sample and the corresponding 
distributions using 7+jets events. 
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Figure 3: Discovery potential in the inclusive 3 jets + JfiT 
final state at CMS. 



3.2. Exclusive two jets channel 

It expected that in the early days the Bj observable 
will be poorly understood. As studied in [7|], a new dis- 
criminating variable (or) is used instead of I^t- the 
ratio between the transverse energy of the second lead- 
ing jet and the transverse invariant mass of the two 
leading jets. The search is performed in events with 
the two leading jets having px> 50 GeV/c, and no 
hard third jet (pr< 50 GeV/c) or hard leptons (pt< 
10 GeV/c). The dominant QCD dijet background can 
be dramatically reduced by requiring ax > 0.55 as it 
can be seen in Fig. 2] An additional cut on the HT > 
500 GeV keeps the SM backgrounds under control. 

The remaining backgrounds as listed in Fig. [4] are 
estimated with a data-driven method using an ABCD- 
like technique where the variables of interest are ax 
and |?/| of the leading jet. As it can be seen in Fig. /ref- 
fig:alphaTeta, the signal (LM1) is expected to have the 
leading jet mostly in the central region (\rj\ < 2.5), 
while the SM background extends more forward in 
pseudo-rapidity. 
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Figure 4: Distribution 
events with twojets. 
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Figure 6: The ratio of the number of events with ar > 
0.55 to the number of events with ot < 0.55 as a function 
of the pseudo-rapidity of the leading jet. The triangles 
correspond to the SM background events while the squares 
represent the mix of signal (LM1) and the various SM 
backgrounds. 




Figure 5: Distribution of the pseudo-rapidity of the lead- 
ing jet for SUSY signal (LM1) and various SM back- 
grounds. 



The signal region (C) is denned by qj- > 0.55 and 
771 1 <2.5, while the control regions are defined as fol- 
lows: region A has ar < 0.55 and >2.5, region 
B has olt > 0.55 and \r]i\ >2.5, and region D has 
cxt < 0.55 and |?7i| <2.5 The amount of background in 
the signal region is then equal to DxB/A. Fig.[6]shows 
that for background (triangles) the fraction B/A is 
constant as a function of the pseudo-rapidity of the 
leading jet, while adding signal (squares) introduces 
a dependence. It has to be noted here that relaxing 
the HT cut results in a larger signal contamination in 
regions A and B which will create larger systematic 
uncertainties on the background estimate. 

With this background estimation method, the sig- 
nal to background ratio is expected to be of order 6 
and SUSY signal to be «44 events for 100 pb _1 of 



data. 

3.2.1 . Generalization to n-jets 

The techniques used in the search performed on the 
sample of events with exactly two jets can be applied 
on the sample with more jets by generalizing the defi- 
nition of olt. This can be done by combining the jets 
into two groups such that the difference (AH™ n ) be- 
tween the HT calculated in each group is minimized. 
The generalized a T is defined as (H T -AH™ m )/2M T , 
where HT is computed using all the jets used in the 
analysis and is their transverse invariant mass. 
In Fig. [7] it is shown that the dominant QCD back- 
ground can be effectively reduced by cutting on this 
generalized variable. 

As detailed in 0] , this analysis has been performed 
in sample of events with 3, 4, 5, and 6 jets in the final 
state. The trigger used to select these events requires 
at LI that the leading jet has Et> 70 GeV and at 
HLT E T > 110 GeV. The jets used in this analysis 
are required to have Et> 50 GeV and |ry| <3, while 
the leading two jets must have Et> 100 GeV and 
the leading jet |?/| <2. In order to eliminate leptonic 
events the muons and electrons in these events should 
have p T < 10 GeV/c. 

The remaining backgrounds arc further reduced by 
requiring the events to have HT > 350 GeV. The con- 
tribution from the QCD background is diminished by 
requiring that q?t > 0.55. To reduce the QCD con- 
tribution to large ax > values, it is required that the 
ratio of MHT calculated using jets with Et> 50 GeV 
to the MHT using jets with E T > 30 GeV is less than 
1.25. 

As in the case of the exclusive two jets analysis, the 
signal to background ratio is expected to be of order 
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Figure 7: Distribution of the clt variable calculated in 
events with n jets. 



6 and SUSY signal to be «52 events for 100 pb 1 of 
data. 

3.3. Conclusion 

The CMS collaboration has a search for SUSY pro- 
gram covering all possible final states and it is using a 
variety of tools and methods to this end. In particular, 



the searches for SUSY in the fully hadronic channels 
are well advanced and the most sensitive. This sen- 
sitivity is in part due to the data-driven approaches 
to estimate the SM backgrounds thus minimizing the 
systematic uncertainties. The studies described in this 
report indicate that SUSY can be discovered in the 
fully hadronic final state with less than 100 pb -1 of 
understood data. 
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